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Abstract 

The diffractive production of vector bosons at the Tevatron is studied. 
We take a hard pomeron flux and use the HI parton density functions 
extracted from their measurement of F®^. To this we add a reggeon ex- 
change contribution. We find that the ratio of diffractive to non-diffractive 
W boson production is in good agreement with the CDF data. We note 
that the poorly understood reggeon contribution might be as large as the 
pomeron contribution in the kinematic range of the data. We also note 
that our pomeron exchange contribution is much smaller than earlier pre- 
dictions due to our use of a hard pomeron flux. All our results have been 
obtained using the lowest order matrix elements and we have compared 
to the results obtained using a modified version of Herwig which is ca- 
pable of generating a wide variety of diffractive scattering processes. Gap 
survival is estimated using Pythia and is found to be around 60%. We 
make predictions which could be compared to future measurements of 
diffractive Drell-Yan production. 
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1 Introduction 



So far there is no consensus on the interpretation of the hard diffractive scat- 
tering processes measured at the Tevatron [Q, §, |||, @]- Invariably the rates 
measured at the Tevatron are much smaller than would be expected using 
the diffractive parton density functions measured in diffractive deep inelastic 
scattering at Hera j^, ||. On the one hand this is not in contradiction with 
theory since the diffractive parton density functions extracted in deep inelastic 
scattering are not expected to be applicable in diffractive scattering involving 
incoming hadrons 0. However, there is the hope that this breaking of univer- 
sality might be described in terms of an overall multiplicative "gap survival" 
factor f|. 

In this paper we turn our attention to diffractive vector boson production 
at the Tevatron. We perform an analysis using the diffractive parton density 
functions extracted by the HI collaboration || with one important modification. 
Namely, we use a larger effective pomeron intercept, i.e. ap(0) = 1.45 rather 
than the value 1.203 extracted by HlQ We believe this to be a natural choice 
given that the scale associated with diffractive W production is some two orders 
of magnitude larger than that in diffractive DIS at Hera, and given that there 
is evidence for an effective pomeron intercept which increases slowly as the 
hard scale increases |], [nj. We include the rather poorly constrained reggeon 
exchange contribution, also as determined by the HI collaboration 

In Section ^ we focus on W boson production as measured by CDF [01. 
This process has been studied within a similar framework previously [ fl2| , [L3| 
where it was found that the rates based upon direct extrapolation of the deep 
inelastic parton density functions typically overestimate the CDF data by an 
order of magnitude. However, neither of these analyses used a hard pomeron 
flux, nor did they consider the reggeon contribution. Possible explanations for 
the largeness of the theoretical cross-section have been presented ||l3] , |T^| , some 
based upon the concept of gap survival Jl5| , fl6l| . We shall discuss the issue of 
gap survival in Section [|. 

Throughout our study, we shall perform our calculations using the relevant 
leading order electroweak matrix elements. In addition, we compare to results 
obtained using a modified version of Herwig, which we call Pomwig, that is 
able to generate diffractive scattering in a wide range of processes jl^, |l9|J^| 

Our studies of diffractive W production conclude that the reggeon contri- 
bution is probably not negligible. It is poorly understood and more studies are 
needed in order to reduce the uncertainty from this source. Furthermore, we 
propose that the rate for pomeron initiated processes is much smaller than has 
been found hitherto, due mainly to our use of diffractive parton densities with 
a steeper Xjp dependence than occurs in softer pomeron exchange models. In 

1 We speak of the "effective pomeron intercept" since it is commonplace. However, what we 
refer to is the parameter which controls the x w distribution, assuming a factorisation of the 
diffractive parton density into a product of a function of x w and a function of (3. 

2 Thp mrtp can hp nhtainpd from the authors upon request or from 
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Section [| we examine diffractive Drell-Yan production in order to investigate 
the additional information that can be obtained at the Tevatron. In Section 
|B] we present our conclusions. 



2 Diffractive W production 

For diffractive vector boson production in proton-antiproton collisions with a 
fast forward proton we use the general formula 

% = TsX< I d{cose) I IT I M F */p( x *)fiMf)fi*( x >f) s ^ W 



for the rapidity (77) distribution of the lepton which is produced in the decay 
of the boson. Using the language of Ingelman & Schlein |20), the pomeron 
flux factor is Fjp/pix^,) and /i,ip(/3, /J, 2 ) is the parton distribution function for 
partons of flavour i in the pomeron. As is conventional, x-^ is the fraction of 
the proton energy carried by the pomeron and f3 is the fraction of the pomeron 
energy carried by the parton involved in the hard scatter. The pp centre-of-mass 
energy is denoted s = (1800 GeV) 2 . fj y p(x,fj, 2 ) is the anti-proton distribution 
for partons of flavour j. Throughout we take the factorisation scale fi to be the 
subprocess centre-of-mass energy s. 

In the particular case of W~ production we use the narrow width approxi- 
mation to write the rapidity distribution of the scattered electror^] (throughout 
this paper positive rapidity is defined in the proton direction): 

% = / d ( cos0 ) J ^^ip/ p (^)/wp(/3,M 2 )/^^M 2 ) Va Q% 

(2) 

Vij = cos 2 C for ud and cs incoming partons or Vij — sin 2 9 C for us and cd 
incoming partons. Qy = u in the case that a quark from the pomeron partici- 
pates in the hard scatter and Qfj = t when it is an anti-quark, where 

u = -^(l + cos6), 
M 2 

£=-^(l-OOBtf), 

and is the scattering angle of the scattered electron in the parton-parton 
centre-of-mass frame measured relative to the proton direction. In the case of 
W + production it is necessary to make the interchange t <-» u if we are to 
define 77 as the rapidity of the scattered positron. We integrate over | cos#| < 
(1 - 4p 2 min /Af^) 1 / 2 with p tmin = 20 GeV and all x w < 0.1. These are the cuts 
used by the CDF collaboration. 

3 For definiteness we speak of the 1st generation of leptons only. Including the 2nd gener- 
ation will of course increase all our cross-sections by a factor of two. 
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In order to compute the ratio of diffractive to non-diffractive events we use 
a similar expression for the non-diffractive cross-section, i.e. by replacing the 
pomeron flux factor by 6(1 — x w ) and the pomeron parton densities by the 
proton ones. In computing this ratio we integrate over the CDF range, \r/\ < 1.1. 

Our default is to use the HI pomeron flux, with an increased value for the 
pomeron intercept, i.e. 

f / 1 \ 1 - 90+0 - 52t 

Fjp /p (x„) = N J dt exp(4.6i) l-\ (3) 

and integrate over all t (in practise it is sufficient to integrate over < — t < 1 
GeV 2 ). The normalisation N is defined so that the flux is equal to 1 at x w = 
0.0030 Note that this is a significantly steeper flux than that proposed by 
Donnachie & Landshoff (dl)^2|, which is roughly equivalent to the above flux 
except for the replacement of 1.90 by 1.17. For the pomeron parton densities we 
use the favoured leading order HI fit 2. We always use the GRV94 LO parton 
densities for incoming hadrons p3[ . 

For the reggeon we simply replace the pomeron flux with |24| 

„ / 1 \ °- 50+0 - 90t 

Fm/pM = ^Jdt exp(4.0i) ( — J (4) 

and the pomeron parton densities are replaced by the GRV leading order pion 
densities of J2^]. This is the reggeon contribution that was shown to agree 
with the HI data collected in the range 0.1 < x w < 0.3 Jn|. We ignore any 
possible interference between the pomeron and the reggeon although this may 
be important. 

In Figure [j] we show the rapidity distribution of the scattered lepton for both 
the pomeron and reggeon contributions. In all cases the electron distribution 
is more forward peaked than the positron due to the chiral nature of the W 
interactions. Note the significant suppression obtained when using the hard flux 
as compared to the DL flux. This is a consequence of the flatter x w distribution 
in the latter case - a small effect at x w values typical for Hera but not so for the 
much larger x w values probed at the Tevatron. The fact that the Tevatron 
is at significantly larger x w is also responsible for the large contribution from the 
reggeon. Also significant is the fact that the peaks of the lepton distributions 
are more backward when using the hard flux compared to the DL flux. This 
arises because the larger x m values in the DL case lead to a smaller boost to 
negative rapidity than occurs using the HI flux. The implication is that when 
we integrate over the CDF acceptance region (\t]\ < 1.1) a larger fraction of 
the leptons produced using the hard flux fall outside the region of acceptance 
thereby still further reducing the predicted cross-section. 

The bottom line is that the ratio of diffractive to non-diffractive events, 
summed over the electron and positron, and multiplied by 2 to account for 
proton and antiproton dissociation, is much smaller in the case of the hard flux 

Using the routine QCD 1994|2l). 
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Figure 1: The rapidity distribution of the electron and positron for the pomeron 
(upper plot) and reggeon (lower plot) contributions in diffractive W production. 
In all cases the electron distribution is peaked at more positive rapidity than 
the positron. In the pomeron case, the broken lines correspond to using the flux 
of Donnachie & Landshoff rather than our default flux. The histogram is the 
result from Pomwig. 

compared to the DL flux. In particular we find i?^ L = 6.38% and R^ ald = 0.86% 
where i? 13 denotes the ratio of pomeron exchange events to non-diffractive events 
integrated over the CDF acceptance region. This is to be compared to the CDF 
value for the ratio of diffractive to non-diffractive events, -Rcdf = 1.15 ± 0.55%. 
However, we cannot avoid the reggeon contribution. Including it leads to a ratio 

iC? = 1-67%. 

Using Pomwig, this ratio is only very slightly different at 1.72%. 

Before leaving this section it is necessary to make a few remarks regarding 
the uncertainty in this theoretical result. There is a large uncertainty associated 
with what we have called the effective pomeron intercept, since there is no 
existing measurement of this parameter in the kinematic range relevant to the 
Tevatron. We have used the HI fits for the diffractive parton density functions. 
These density functions are only evolved up to Q 2 = 75 GeV 2 and remain 
constant thereafter. We have also obtained results using the ACTW parton 
density functions fl2]| , which are evolved to high Q 2 and from which we are 
able to estimate that normal QCD evolution will change our prediction for the 
gap fraction by no more than 20%. Next-to- leading order QCD corrections are 
expected to play a role at the 10-20% level. There is also an uncertainty in the 
size of the reggeon contribution that is hard to quantify. 
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3 Gap survival 



So far we have ignored the fact that factorisation is not expected to hold in hard 
diffraction at the Tevatron. One might attempt to account for this physics by 
an overall multiplicative factor which specifies the probability that the rapidity 
gap is destroyed by additional interactions which are not strongly correlated 
with the primary hard scatter. We aim in this section to add plausibility to this 
assumption and to make an estimate for the gap survival factor. 

In a previous paper |^6| we used the multiple interactions model implemented 
in the Pythia event generator to estimate the gap survival probability in gap- 
between jets events |2S|] . In this model the number of additional scatters, 
given a particular hard trigger process, is dependent only on the total energy of 
the collisions and of the scale of the hard trigger process. The estimate of the 
survival probability obtained was 22% for 1800 GeV and 35% for 630 GeV and 
was largely independent of the transverse momentum and rapidity of the jets. 

The increase of multiple interactions with the total energy is simply due to 
the small-x proliferation of gluons which are available for additional scatterings. 
The dependence on the hard scale can be described in a simplified manner as 
follows. 

The model assumes that the partons in the proton are distributed accord- 
ing to a double Gaussian in impact-parameter space. This means that central 
collisions with a small impact parameter have a higher probability for multiple 
scatterers than do peripheral ones with a large impact parameter. For the case of 
VF-production, which we are interested in here, the procedure is to generate the 
hard electroweak scattering and then to choose an impact parameter depending 
on the scale of the process - the higher the scale the larger the probability to 
pick a small impact parameter. The overlap is then calculated from this and the 
number of additional scatterers is obtained. (Note that the procedure becomes 
more complicated if a trigger process with large cross section is used, since the 
eikonalisation of the cross section becomes important in that case.) 

In principle it is possible to run Pompyt fllT] with the multiple interactions 
of Pythia. This will, however, only give multiple interactions between the 
colliding anti-proton and pomeron, and such scatterings will not destroy the gap. 
What we need is to estimate the multiple interactions between the proton and 
anti-proton. Here the scale which determines the overlap is not the hard scale 
of the electroweak lU-production process, but rather the momentum transfer of 
the pomeron — t, which is small. So, although the overlap between the pomeron 
and anti-proton is large, the overlap between the proton and anti-proton may be 
small, resulting in a gap-survival probability much larger than the ones quoted 
above. 

Since it is not possible simply to add multiple interactions between the pro- 
ton and anti-proton in Pompyt, we use the fact that the amount of multiple 
scattering only depends upon the scale of the process and the total collision 
energy. Using different triggering processes which are able to go to small scales 
without large eikonalisation effects in Pythia (using default parameters) we 
find that the probability not to have any additional scatters saturates for small 
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Figure 2: The fraction of events without additional scatterings (Pq) as a func- 
tion of the hard scale in photon exchange (full lines) and di-photon production 
(dashed lines) events. The upper curves are the fraction of events with no ad- 
ditional colourful scatterings whilst the lower ones are the fractions with no 
additional scatterings at all. 

momentum transfers at around 50%. In Figure § we show this probability for 
proton-antiproton scattering via the electromagnetic qq' — > qq' and qq — » 77 
sub-processes. 

For large scales the probability is well below the 22% we quoted in J2(|. This 
is because of the way partons from multiple scatterings are colour-connected to 
the rest of the event in Pythia. With the default parameters only 1/3 of the 
multiple scatterings will be colour-connected to the rest of the event, while the 
remaining scatterings will produce colour singlet qq or gg-pairs which will not 
necessarily destroy the gap. In Figure || we therefore also show the probability 
of having at least one colourful secondary scatter. The true gap fraction ought 
to lie somewhere in between these two extremes. 

There are large uncertainties here. Multiple interactions are poorly under- 
stood and the multiple interaction model in Pythia has not been tuned to 
describe these processes. It is, however, natural to assume that any rapidity 
gap produced due to the exchange of a pomeron or a reggeon may be destroyed 
by multiple interactions. It is also very intuitive that these processes are more 
peripheral than other hard interactions and that the survival probability there- 
fore should be larger than in the case of, e.g. high-i diffraction. The model 
in Pythia does attempt to take such physics into account. Nevertheless, it is 
clear that more effort must be devoted to understanding these processes before 
we can start to compare the diffractive results from Hera with the ones from 
the Tevatron with confidence. 

The results of this section suggest that we should take a gap survival factor 
of around 60% for low t diffraction at the Tevatron, running at 1.8 TeV. This 
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Figure 3: The cross-section as a function of invariant mass for the Drell-Yan 
process. The solid line is the pomeron contribution and the dashed line is the 
reggeon contribution. The curves are obtained for pr > 10 GeV. The dotted line 
is the pomeron contribution for a smaller effective intercept and the dot-dash 
line shows the same but with a variable effective intercept, as explained in the 
text. The histrograms are the Pomwig result. 



should be compared to other models which typically predict gap survival factors 
of 10% or less |fl^ 16 1. Accounting for this gap survival factor we find a gap 
fraction of 



R 



P+]R 
hard 



1.00% 



which is in good agreement with the Tevatron measurement. 



4 Diffractive Drell-Yan production 

In order to test our understanding of diffractive vector particle production, it 
would be useful not only to measure the rapidity distribution of the charged 
lepton in production but also to do the same for the diffractive Drell-Yan 
production process and in this section we present our predictions. 
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The rapidity distribution of the scattered lepton can be computed: 



The advantage of the Drell-Yan process is the ability to tune the invariant mass 
of the dilepton pair and in Figure ||| we show the cross-section as a function of 
invariant mass, M = y/I, integrated over —1.1 < 77 < 1.1 for both the pomeron 
and reggeon contributions. The solid black line represents a hard pomeron with 
effective intercept aip(O) = 1.45, the dotted line represents a hard pomeron 
with an effective intercept ajp(O) = 1.20, and the broken line is the reggeon 
contribution. The histrograms are the Pomwig results, and the non-vanishing 
at low M is a result of the correct matching to the NLO matrix elements that is 
performed in Herwig 6.1. 

To match smoothly with the Hera data it is perhaps more reasonable to 
use an effective intercept that depends upon the hard subprocess scale. This 
possibility leads to the dash-dot curve in Figure || where we have used 

oip(O) = 1.10 + 0.038 InM 2 . 

We wish to stress that this is an entirely phenomenological parameterisation. 
From a theoretical perspective, such an x w dependence cannot be generated by 
exchange of a single regge pole. However this effective behaviour could describe 
the interplay of two or more regge exchanges. For example, in a two-pomeron 
model one might imagine that the harder pomeron, with intercept around 1.45, 
dominates at the Tevatron whilst at Hera there is an important admixture 
of the softer pomeron - this would lead to qualitatively the behaviour we ex- 
amine here |l(J. Alternatively, within the leading-twist PQCD formalism, the 
x w -dependence should not exhibit any Q 2 dependence since all of the evolution 
is contained within what we call the pomeron parton distribution functions. In 
this case, one should interpret our parameterization as an indication that there 
is an important higher-twist component in the relatively low Q 2 Hera data 
that is absent in the Tevatron data. 



5 Conclusions & Outlook 

We have shown that the current Tevatron data on diffractive W ± production 
can be understood in terms of the diffractive parton densities extracted from 
Hera data provided we use a hard pomeron intercept, ap(0) = 1.45. Modelling 
the diffractive parton densities in terms of pomeron and reggeon contributions, 
we find that the reggeon contribution is approximately equal to the pomeron 
contribution in the kinematic regime relevant for the Tevatron. Moreover, we 
propose that the pomeron contribution is much smaller than previous estimates. 
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We used Pythia to model the gap survival probability and found it to be quite 
large at 60%. In conclusion we find a gap fraction in production of 1.00% 
which is to be compared to the CDF result of (1.15 ± 0.55)%. 

A precise comparison between theory and data is not yet possible due pri- 
marily to lack of understanding of gap survival effects and of the reggeon contri- 
bution. Future data collected on diffractive Drell-Yan production at the Teva- 
TRON will certainly help clarify the physics, not least by providing valuable 
data over a wide range of scales, from those typical of dis at Hera to those 
in diffractive Z° production. Measurements to lower x w at the Tevatron will 
also be important in helping disentangle the competing effects. 

The natural next step is to compare the model presented here to the data 
for diffractive dijet production [|l| and diffractive central jet production ||. Pre- 
liminary investigations suggest that the existing contradictions will be resolved 
- we deal with these processes in a future publication. 

Acknowledgements 

We should like to thank Dino Goulianos, Hannes Jung, Andy Mehta, Paul 
Newman, Julian Phillips and Mike Seymour for their valuable assistance. 

References 

[1] F. Abe et al., CDF Collaboration, Phys. Rev. Lett. 79, 2636 (1997); 
T. Affolder et al., CDF Collaboration, Phys. Rev. Lett. 84, 5043 (2000) 

[2] F. Abe et al., CDF Collaboration, Phys. Rev. Lett. 78 (1997) 2698. 

[3] T. Affolder et al., CDF Collaboration, Phys. Rev. Lett. 85, 4215 (2000). 

[4] T. Affolder et al., CDF Collaboration, Phys. Rev. Lett. 84, 232 (2000). 

[5] T. Ahmed et al., HI Collaboration, Phys. Lett. B348 (1995) 681; 
C. Adloff et al., HI Collaboration, Z. Phys. C76 (1997) 613. 

[6] J. Breitweg et al., ZEUS Collaboration, Eur. Phys. J. CI (1998) 81; 
Eur. Phys. J. C6 (1999) 43. 

[7] J.C. Collins, Phys. Rev. D57 (1998) 3051. 

[8] J.D. Bjorken, Phys. Rev. D47 (1993) 10. 

[9] For example see B.E. Cox, J. Phys. G25 (1999) 1377. 

[10] A. Donnachie and P.V. Landshoff, Phys. Lett. B437 (1998) 408. 

[11] C. Adloff et al., HI Collaboration, Eur. Phys. J. C6 (1999) 587. 

[12] L. Alvero, et al., Phys. Rev. D59 (1999) 074022 



10 



[13] R.J.M. Covolan and M.S. Soares, Phys. Rev. D60 (1999) 054005; 
erratum-ibid. D61 (2000) 01990. 

[14] K. Goulianos, Phys. Lett. B358 (1995) 379; 
C.I. Tan, Phys. Rep. 315 (1999) 175. 

[15] E. Gotsman, E. Levin and U. Maor, Phys. Rev. D60 (1999) 094011. 

[16] V.A. Khoze, A.D. Martin and M.G. Ryskin, 

\protect\vrule widthOpt\protect\href {http : //arXiv . org/abs/hep-ph/0007083}{hep-ph/000708: 

[17] P. Bruni and G. Ingelman, DESY 93-187 (1993). 
www3 . tsl .uu. se/thep/pompyt/ 

[18] G. Marchesini et al., Comput. Phys. Commun. 67 (1992) 465; 

\protect\vrule widthOpt\protect\href {http : //arXiv . org/abs/hep-ph/0011363Mhep-ph/001136; 

[19] B.E. Cox and J.R. Forshaw, \pr otect\vrule widthOpt\protect\href {http : / /arXiv . org/abs/hep-ph, 
www . hep . man . ac . uk/ pomwig/ 

[20] G. Ingelman and P. Schlein, Phys. Lett. B152 (1985) 256. 

[21] H. Jung and J. P. Phillips, private communication. 

[22] A. Donnachie and P.V. Landshoff, Phys. Lett. B191 (1987) 309; 
Nucl. Phys. B303 (1988) 634. 

[23] M. Ghick, E. Reya and A. Vogt, Z. Phys. C67 (1995) 433. 

[24] K. Golec-Biernat, et al., Phys. Rev. D56 (1997) 3955; 

A. Szczurek, N.N. Nikolaev and J. Speth, Phys. Lett. B428 (1998) 383. 

[25] M. Gliick, E. Reya and A. Vogt, Z. Phys. C53 (1992) 651. 

[26] B.E. Cox, J.R. Forshaw and L. Lonnblad, JHEP (1999) 9910. 

[27] T. Sjostrand, Comput. Phys. Commun., 82 (1994) 74. 
www. thep . lu. se/~torbjorn/Pythia.html 

[28] T. Sjostrand and M. van Zijl, Phys. Rev. 36 (1987) 2019 



11 



